VARIATIONAL APPROACH TO COMPLICATED SIMILARITY 
SOLUTIONS OF HIGHER-ORDER NONLINEAR PDES. I 

V.A. GALAKTIONOV, E. MITIDIERI, AND S.I. POHOZAEV 

Abstract. The Cauchy problem for (x,t) E M.^ x M-|_ for three higher-order degenerate 
quasiUncar PDEs, as basic models, 

ut = (-1)'"+iA™(|m|"m) + 
utt = (-l)™+iA'"(|u|"u) + 
Ut = (-l)™+MA"(|u|"u)].i + {\u\"u),„ 

where n > is a fixed exponent and A™ is the (m > 2)th iteration of the Laplacian, 
is studied. This diverse class of degenerate PDEs embraces equations of different three 
types: parabolic, hyperbolic, and nonlinear dispersion. Such degenerate evolution equa- 
tions from various applications in mechanics and physics admitting compactly supported 
and blow-up solutions attracted attention of the mathematicians since the 1970-80s. 

In the present paper, some general local, global, and blow-up features of such PDEs 
on the basis of construction of their blow-up similarity and travelling wave solutions are 
revealed. Blow-up, i.e., nonexistence of global in time solutions is proved by various 
methods. In particular, for m = 2 and m = 3, such similarity patterns lead to the 
following semilinear fourth- and sixth-order elliptic PDEs with non-coercive operators 
and non-Lipschitz nonlinearities: 

(0.1) -A^i^-I-F- |F|"^F = and A^i^ -|- F - in x M+, 

which were not addressed before in the mathematical literature. The goal is, using a 
variety of analytic variational, qualitative, and, often, numerical methods, to justify that 
equations (|0.ip admit an infinite at least countable set of countable families of compactly 
supported solutions that are oscillatory near finite interfaces. In a whole, this solution 
set exhibit typical features of being of a chaotic structure. 

The present paper is an earlier extended version of [17] . In particular, here we pay 
more attention to some category/fibering aspects of critical values and points, as well as 
to sixth and higher-order equations with m > 3, while in [T7] the case m — 2 was under 
maximally detailed scrutiny. 



1. INTRODUCTION: HIGHER-ORDER MODELS AND BLOW-UP/COMPACTON SOLUTIONS 

1.1. Three types of nonlinear PDEs under consideration. We describe common 
local and global properties of weak compactly solutions of classes of nonlinear partial 
differential equations (PDEs) of parabolic, hyperbolic, and nonlinear dispersion type: 
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(I) 2mth-order quasilinear parabolic equations with regional blow-up, 

(II) 2mth-order quasilinear hyperbolic equations with regional blow-up, and 

(III) (2m+l)th-order nonlinear dispersion equations with compactons. 

Therefore, we plan to study some common local and global properties of weak com- 
pactly solutions of three classes of quasilinear partial differential equations (PDEs) of 
parabolic, hyperbolic, and nonlinear dispersion type, which, in general, look like having 
nothing in common. Studying and better understanding of nonlinear degenerate PDEs 
of higher-order including a new class of less developed nonlinear dispersion equations 
(NDEs) from compacton theory are striking features of modern general PDE theory at 
the beginning of the twenty first century. It is worth noting and realizing that several key 
theoretical demands of modern mathematics are already associated and connected with 
some common local and global features of nonlinear evolution PDEs of different types and 
orders, including higher-order parabolic, hyperbolic, nonlinear dispersion, and others as 
typical representatives. 

Regardless the great progress of PDE theory achieved in the twentieth century for many 
key classes of nonlinear equations [B], the transition process to higher-order degenerate 
PDEs with more and more complicated non-monotone, non-potential, and non-symmetric 
nonlinear operators will require different new methods and mathematical approaches. 
Moreover, it seems that, for some types of such nonlinear higher-order problems, the 
entirely rigorous "exhaustive" goal of developing a complete description of solutions, their 
properties, functional settings of problems of interest, etc., cannot be achieved in principle, 
in view of an extremal variety of singular, bifurcation, and branching phenomena that are 
contained in such multi-dimensional evolution. In many cases, the main results should 
be extracted by a combination of methods of various analytic, qualitative, and numerical 
origins. In many cases, the main results should be extracted by a combination of methods 
of various analytic, qualitative, and numerical origin^. 

In the present paper, we deal with complicated pattern sets, where, for the elliptic 
problems in and even for the corresponding one-dimensional ODE reductions, using 
the proposed analytic-numerical approaches is necessary and unavoidable. As a first 
illustration of such features, let us mention that, according to our current experience, for 
such classes of second-order variational problems, 

distinguishing the classic Lusternik-Schnirel'man countable sequence of 
critical values and points is not possible without refined numerical methods, 

in view of huge complicated multiplicity of other admitted solutions. It is essential, that 
the arising problems do not admit, as customary for other classes of elliptic equations, 
any homotopy classification of solutions (say, on the hodograph plane), since all the com- 
pactly supported solutions are infinitely oscillatory that makes the homotopy rotational 
parameter infinite and hence the method non- applicable. 

""^This is not a novelty in modern mathematics, where several fundamental rigorous results have been 
already justified with the aid of hard, refined, and reliable numerical experiments; q.v. e.g., Tucker's 
proof of existence of a robust strange attractor for the 3D Lorenz system [ID] . 
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Let us now introduce these three classes, (I)— (III), of PDEs and corresponding non- 
linear phenomena to be studied by some unified approaches. 



1.2. (I) Combustion models: regional blow-up, global stability, main goals, 
and first discussion. We begin with the following quasilinear degenerate 2mth-order 
parabolic equation of reaction-diffusion (combustion) type: 

(1.1) ut = (-l)™+^A™(|ii|"n) + in X M+, 

where n > is a fixed exponent, m > 2 is integer, and A denoted the Laplace operator 
in R. 

Globally asymptotically stable exact blow-up solutions of S-regime. In the simplest case 
m = 1 and = 1, ( 11. ip is nowadays the canonical quasilinear heat equation 

(1.2) Ut = {u'^-^^)^^ + in M X M+ {u > 0), 

which occurs in combustion theory. The reaction-diffusion equation (11. 2p . playing a key 
role in blow-up PDE theory, was under scrutiny since the middle of the 1970s. In 1976, 
Kurdyumov, with his former PhD students, Mikhailov and Zmitrenko [q.v. [3H] and [371 
Ch. 4] for history) discovered the phenomenon of heat and combustion localization by 
studying the blow-up separate variables Zmitrenko-Kurdyumov solution of (11.21) : 

(1.3) us{x,t) = {T -ty^ f{x) in Mx(0,r), 
where T > is the blow-up time, and / satisfies the ODE 

(1.4) i / = (r+i)" + for xeR. 

It turned out that (11.41) possesses the exphcit compactly supported solution 



[St§cos^(i(^)]^ if N<^ 



'1.5) fix) = < Ln(n+2) l2(n+l). 

This explicit integration of the ODE (ll.4p was amazing and rather surprising in the middle 
of the 1970s and led then to the foundation of blow-up and heat localization theory. In 
dimension > 1, the blow-up solution (11.31) does indeed exist [23 P- 183] but not in an 
explicit form (so that, it seems. [T31) is the only available elegant form). 

Blow-up S-regime for higher-order parabolic PDEs. Evidently, the 2mth-order counter- 
part (11.11) admits the regional blow-up solution of the same form (11. 3p . but the profile 
/ = /(?/) then solves a more complicated ODE 

(1.6) (-l)™+iA-(|/r/) + l/r/ = ^ / in M^. 

After natural change, this gives the following equation with a non-Lipschitz nonlinearity: 

F=\f\^f =^ (-l)™+iA'"F + F- i |F| "+iF = in M^. 
Finally, we scale out the multiplier - in the nonlinear term. 



;i.7) F^n-'^F 



-lyn+l^rap p n+l _p = Q 



m 



In the one-dimensional case = 1, we obtain a simpler ODE, 



(1.8) 



F n 



n+l 



F 




+^F = in M. 



Thus, according to (11. 3p . the elliptic problems ( II. 7p and the ODE ( II. 8p for = 1 are 
responsible for the possible "geometrical shapes" of regional blow-up described by the 
higher-order combustion model (11.11) . 

Plan and main goals of the paper related to parabolic PDEs. Unlike the second-order case 
(II. 5p . exp/zcit compactly supported solutions F{x) of (II. 8p for any m>2 are not available. 
Moreover, it turns out that such profiles F{x) have rather complicated local and global 
structure. We are not aware of any rigorous or even formal qualitative results concerning 
existence, multiplicity, or global structure of solutions of ODEs such as (11.81) . Our main 
goals are four-fold: 

(ii) Problem "Blow-up" : proving finite-time blow-up in the parabolic (and hyper- 
bolic) PDEs under consideration (Section [2]); 

(ii) Problem "Multiplicity" : existence and multiplicity for elliptic PDEs (11.71) and 
the ODEs ([H]) (Section [3]); 

(iii) Problem "Oscillations" : the generic structure of oscillatory solutions of (11.81) 
near interfaces (Section Hj); and 

(iv) Problem "Numerics" : numerical study of various families of F{x) ( Sections [5l [6]). 

The research will be continued in the second half of this paper [18], where we intend 
to refine our results on the multiplicity of solutions (especially, for m > 2), pose the 
problem on a "Sturm index" of solutions (a homotopy classification of some sub-families 
of solutions), and introduce and study related analytic models with similar families of 
solutions. 

Thus, in particular, we show that ODEs (II. 8p . as well as the PDE ( II. 7p . for any m > 2 
admit infinitely many countable families of compactly supported solutions in M, and 
the whole solution set exhibits certain chaotic properties. Our analysis will be based 
on a combination of analytic (variational and others), numerical, and various formal 
techniques. Explaining existence, multiplicity, and asymptotics for the nonlinear problems 
involved, we leave several open mathematical problems. Some of these for higher-order 
equations are extremely difficult. 

1.3. (II) Regional blow-up in quasilinear hyperbolic equations. Secondly, con- 
sider the 2mth-order hyperbolic counterpart of (II. ip . 



(1.9) Utt = (-l)'"+^A"^(|ii|"ix) + lul^u in x 
We begin the discussion of its blow-up solutions in ID, i.e., for 

(1.10) Utt = (m"+^)xx + u""^^ in M X M+ (u > 0). 
Here the blow-up solutions and the ODE take the form 



(1.11) us{x,t) = iT-t)-^^fix) =^ f(f + i)/ = (r+^) 
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Using extra scaling, 

(1-12) /»= 

yields the same ODE fll.4p and hence the exact localized solution fll.51) . 

For the A^-dimensional PDE fll.9p . looking for the same solution (11.111) . after scaling, 
leads to the elliptic equation (11. 7p . 

1.4. (Ill) Nonlinear dispersion equations and compactons. In a general setting, 
these rather unusual PDEs take the form 

(1.13) Ut = {-ir+^[/y^{\uY'u%, + {\u\''u),, in M^xM+, 

where the right-hand side is just the derivation D^^ of that in the parabolic counterpart 
(II. ip . Then the elliptic problem (11.70 occurs when studying travelling wave (TW) solutions 
of (ll.lSp . As usual, we explain this first in a simpler ID case. 

Let = 1 and m = 1 in (11.130 that yields the third-order Rosenau-Hyman (RH) 
equation 

(1.14) Ut = {U^)xxx + {u^)x, 

which is known to model the effect of nonlinear dispersion in the pattern formation in 
liquid drops [.36j. It is the K{2, 2) equation from the general K{m, n) family of nonlinear 
dispersion equations (NDEs) 

(1.15) = (u")x.. + (Ox («>0), 

that also describe phenomena of compact pattern formation, [331 El] • In addition, such 
PDEs appear in curve motion and shortening flows [35]. Similar to the previous models, 
the K{m, n) equation (ll.lSp with n > 1 is degenerated at -u = 0, and therefore may exhibit 
finite speed of propagation and admit solutions with finite interfaces. A permanent source 
of NDEs is integrable equation theory, e.g. look at the integrable fifth-order Kawamoto 
equation [23j (see [201 Ch. 4] for other models), which is of NDE's type: 

(1.16) 

Questions of local existence, uniqueness, regularity, shock and rarefaction wave formation, 
finite propagation and interfaces, including degenerate higher-order models to be studied 
are treated in [T9l [Ti] : see also comments in [20l Ch. 4.2]. 

We will study some particular continuous solutions of the NDEs that give insight on 
several generic properties of such nonlinear PDEs. The crucial advantage of the RH equa- 
tion (I1.14P is that it possesses explicit moving compactly supported soliton-type solutions, 
called compactons [36], which are travelling wave (TW) solutions: 

Compactons: manifolds of TWs and blow-up S-regime solutions coincide. Let us show that 
such compactons are directly related to the blow-up patterns presented above. Actually, 
explicit TW compactons exist for the nonlinear dispersion KdV-type equations with ar- 
bitrary power nonlinearities (formulae will be given shortly) 

(1.17) = (m"+1),,, + (m"+1), in MxM+. 
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This is the K{1 + n, 1 + n) model, |36j . 

Thus, compactons as solutions of the equation fll.171) have the usual TW structure 



(1.18) ixc(x, t) = f{y), y = x- Xt, 
so that, on substitution, / satisfies the ODE 

(1.19) - xf = ir^T + ir^'r, 

and, on integration once, 

(1.20) _A/ = (r+y + r+' + A 

where D G M is a constant of integration. Setting D = 0, which means the physical 
condition of zero flux at the interfaces, leads to the blow-up ODE (11. 4p . so that the 
compacton equation (11.201) coincide with the blow-up one (ll.4p if 

(1.21) -A = i (or -A = f (f + 1) to match dm])). 

This yields the compacton solution ( ll.lSp with the same compactly supported proflle (11. Sp 
with the translation x ^ y = x — Xt. 

Therefore, in ID, the blow-up solutions (II. 3p . (1 1.1 ID of the parabolic and hyperbolic 
PDEs and the compacton solution (11.180 of the nonlinear dispersion equation (11.170 are 
essentially of a similar mathematical (both the ODE and PDE) nature, and, possibly, more 
than that. This reflects a certain universality principle of compact structure formation 
in nonlinear evolution PDEs. Stability features of the TW compacton (I1.18P in the PDE 
setting (I1.17P are unknown, as well as for the higher-order counterparts to be posed next. 

In the A^- dimensional geometry, i.e., for the PDE (11.130 . looking for a TW moving in 
the Xi-direction only, 

(1.22) Mc(x,t) = /(Xi - At,X2, ...,XAr) [X = -\) 

we obtain on integration in xi the elliptic problem (11.70 . 

Thus, we have introduced the necessary three classes, (I), (II), (HI), of nonlinear 
higher-order PDEs in X M , , which, being representatives of very different three equa- 
tions types, nevertheless will be shown to exhibit quite similar evolution features (if neces- 
sary, up to replacing blow-up by travelling wave moving), and the coinciding complicated 
countable sets of evolution patterns. These common features reveal an exiting concept of 
a certain unifled principle of singularity formation phenomena in general nonlinear PDE 
theory, that we seem just begin to touch and study in the twenty flrst century. Sev- 
eral classic mathematical concepts and techniques successfully developed in the twentieth 
century including, of course, Sobolev legacy continue to be key, but also new ideas from 
different ranges of various rigorous and qualitative natures are desperately needed for 
tackling such fundamental difficulties and open problem arising. 
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2. Problem "Blow-up": GENERAL BLOW-UP ANALYSIS OF PARABOLIC AND 

HYPERBOLIC PDES 

2.1. On global existence and blow-up in higher-order parabolic equations. We 

begin with the parabolic model fll.ll) . Bearing in mind a compactly supported nature of 
the solutions under consideration, we consider (11 .ip in a bounded domain Q C with 
a smooth boundary dQ, with Dirichlet boundary conditions 

(2.1) u = Du = ... = D'^'-^u = on dVL x M+, 
and a given sufficiently smooth and bounded initial data 

(2.2) u{x,{)) = uq{x) in VL. 

We will show that the phenomenon of blow-up depends essentially on the size of the 
domain. Beforehand, let us observe that the diffusion operator on the right-hand side 
in (II. ip is a monotone operator in H~'^{Q), so that the unique local solvability of the 
problem in suitable Sobolev spaces is covered by classic theory of monotone operators; 
see Lions' book [221 Ch. 2]. We next show that, under certain conditions, some of these 
solutions are global in time but some ones cannot be globally extended and blow-up in 
finite time. 

For convenience, we use the natural substitution 

(2.3) V = \u\^u =^ Vo{x) = \uo{x)\'^uo{x), 

that leads to the following parabolic equation with a standard linear operator on the 
right-hand side: 

(2.4) {ij{v))t = {-l)'^^^A'^v + v in x R+, with ij{v) = \v\~^v, 

where v satisfies the same Dirichlet boundary condition (12.11) . 

Multiplying (12.41) by v in L^(f2) and integrating by parts by using (12.11) yields 

(2.5) ^Tt I \v\^dx = -J \D"'v\'^dx + J v^dx = E{v), 

n n n 

where we use the notation D™ = for even m and -D™ = VA^^ for odd m. By 
Sobolev's embedding theorem, H"^{Q) C L^(fi) compactly, and moreover, the following 
sharp estimates holds: 

(2.6) J v^dx < ^ f \D"'v\'^dx in H^{n), 

where Ai = Ai(fi) > is the ffist simple eigenvalue of the poly-harmonic operator (— A)*" 
with the Dirichlet boundary conditions (12. ip : 

(2.7) (-A)'"ei = Aid in Q, eieH^'^iQ). 

For m = 1, since (—A) > is strictly negative in the metric of L^(fi), we have, by 
Jentzsch's classic theorem (1912) on the positivity of the ffist eigenfunction for linear 
integral operators with positive kernels, that 

(2.8) ei{x) > in Q. 

7 



For m > 2, fl2.8p remains valid, e.g., for the unit ball Q = Bi. Indeed, in the case 
of r2 = Bi, the Green function of the poly-harmonic operator (—A)'" with Dirichlet 
boundary conditions is positive; see first results by Boggio (1901-05) [11,15] (see also Elias 
[10] for later related general results). Again, by Jentzsch's theorem, (12.81) holds. 
It follows from and I^M) that 

^,iJ\vf^dx<{^^-l)J\D-^v\Mx. 
n n 

Global existence for Ai > 1. Thus, we obtain the following inequality: 

(2.9) S ^ / l^^l^dx + (l - ^) / \D"'v\'^dx < 0. 
Consequently, for 

(2.10) \i{n) > 1, 

(12. 6p yields good a priori estimates of solutions in i7 x (0, T) for arbitrarily large T > 0. 
Then, by the standard Galerkin method [26, Ch.l ], we get global existence of solutions 
of the initial-boundary value problem (IB VP) (12.41) . (12.11) . (12. 2p . This means no finite- 
time blow-up for the IBVP provided (I2.10p holds, meaning that the size of domain being 
sufficiently small. 

Global existence for Ai = 1. Note that for Ai = 1, (12. 9p also yields an a priori uniform 
bound, which is weaker, so the proof of global existence becomes more tricky and requires 
extra scaling to complete (this is not directly related to the present discussion, so we omit 
details). In this case, we have the conservation law 

(2.11) f i/j^vi^t)) cidx = Co = Jip{vo)eidx for all t > 0, 

so that by the gradient system property (see below), the global bounded orbit must 
stabilize to a unique stationary solution, which is characterized as follows (recall that Ai 
is always a simple eigenvalue, so the eigenspace is ID): 

(2.12) v{x,t) ^ CqCi^x) as t +oo, where / '0(Coei)ei dx = cq. 

n 

Blow-up for Ai < 1. Let us now show that for the opposite inequality 

(2.13) Al(^]) < 1, 
the solutions blow-up in finite time. 

Blow-up of nonnegative solutions for m = 1. We begin with the simpler case m = 1, 
where, by the Maximum Principle, we can restrict to the class of nonnegative solutions 

(2.14) V = v{x,t) > 0, i.e., assuming that Uq{x) > 0. 

In this case, we can directly study the evolution of the first Fourier coefficient of the 
function ip{v{-,t)). To this end, we multiply (12.40 by the positive eigenfunction ci in 
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L'^{Q) to obtain that 

(2.15) ^ / ^(v)eidx = (1 - Ai) / veidx. 

n n 

In view of (12.141) . we apply Holder's inequality in the right-hand side of (12.151) to derive 
the following ordinary differential inequality for the Fourier coefficient: 

(2.16) ^ > (1 - Ai)c2J"+\ where J{t) = J {x,t)ei{x) dx, ca = ( / eida;)"". 

n n 

Hence, for any nontrivial nonnegative initial data 

Uo{x) ^0 =^ Jo = In vqCi dx > 0, 

we have finite-time blow-up of the solution with the following lower estimate on the Fourier 
coefficient: 

(2.17) J{t)>A{T-t)-v., where A = T = 

On unbounded orbits and blow-up for m >2. It is curious that we do not know a similar 
simple proof of blow-up for the higher-order equations with m > 2. The main technical 
difficulty is that the set of nonnegative solutions (12.141) is not invariant of the parabolic 
flow, so we have to deal with solutions v{x,t) of changing sign. Then, (12.161) cannot be 
derived from (I2.15P by the Holder inequality. Nevertheless, we easily obtain the following 
result first step to blow-up of the orbits: 

Proposition 2.1. Let m>2, {\2.13^ hold, and 

(2.18) E{vo) > 0. 

Then the solution of the IBVP [2^], ^2.1\\ . ^2.2\ is not uniformly bounded for t > 0. 

Proof. We use the obvious fact that ( 12. 4p is a gradient system in H^{fl). Indeed, multi- 
plying (12. 4p by Vt yields, on sufficiently smooth local solutions, 

(2.19) lTtE{v{t)) = T^J\vr^{v,ydx > 0. 

n 

Therefore, under the hypothesis (12.181) we have from (12.51) that 

(2.20) E{v{t)) > E{vo) > =^ ^^fJ\v\'^^dx = E{v)>E{vo)>0, i.e., 

n 

(2.21) J\v{t)\^dx>^E{vo)t^+oo as t ^ +oo. □ 
n 

Concerning the hypothesis (12.181) . recall that by classic dynamical system theory |21j . 
the w-limit set of bounded orbits of gradient systems consists of equilibria only, i.e., 

(2.22) uj{vo) C S = {V e H^{Q) : -{-A)'^V + V = O] . 

Therefore, stabilization to a nontrivial equilibrium is possible if 

A; = 1 for some I >2. 
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Otherwise, we have that 



(2.23) S = {0} (A; ^ 1 for any / > 1). 

Then, formally, by the gradient structure of (12.41) . one should take into account solutions 
that decay to as t — +oo. One can check that (at least formally, a necessary functional 
framework could take some time), the trivial solution has the empty stable manifold, 
so that, under the assumption (12.231) . the result in Proposition 12.11 is naturally expected 
to be true for any nontrivial solution. 

Thus, we have that in the case (12.181) . i.e., for sufficiently large domain Q, solutions 
become arbitrarily large in any suitable metric, including H^{Q) or the uniform one 
Cq{Q). Then, it is a technical matter to show that such large solutions must next blow-up 
in finite time. In fact, often, this is not that straightforward, and omitting this blow-up 
analysis, we would like to attract the attention of the interested Reader. 

Blow-up for m > 2 in a similar modified model. On the other hand, the previous proof of 
blow-up is easily adapted for the following slightly modified equation (12. 4p : 

(2.24) {^iv))t = {-ir^'A"'v + \v\, 

where the source term is replaced by \v\. Actually, for "positively dominant" solutions 
(i.e., for those of a non-zero integral J u{x, t) dx), this is not a big change, and most of our 
self-similar patterns perfectly exist for (12.241) and the oscillatory properties of solutions 
near interfaces remain practically untouched (since the source term plays no role there). 

Take Q = Bi, so that (12.81) holds. Then, instead of (I2.15p . we will get a similar 
inequality, 

(2.25) ^ / ip{v)eidx = J |i;|eidx - Ai / veidx > (1 - Ai) / |-u|eida; > 0, 

n n n n 

where J{t) is defined without the positivity sign restriction, 

(2.26) J{t) = J^{\v\~^v){x,t)ei{x)dx. 
It follows from fl^T^ that, for Ai < 1, 

(2.27) J(0) > ^ J{t) > for t > 0. 
Therefore, by the Holder inequality, 

(2.28) / \v\ei dx > C2{j \v\'^ei dx)"'~^^ > C2(/ \v\~'^vei dx)"'^^ = C2J""^^ 

This allows us to get the inequality (12.161) for the function (I2.26p . Hence, the blow-up 
estimate (I2.17p holds. 



2.2. Blow-up data for higher-order parabolic and hyperbolic PDEs. We have 
seen above that, in general, blow-up occurs for some initial data, since in many cases, 
small data can lead to globally existing sufficiently small solutions (of course, if has a 
nontrivial stable manifold). 
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Below, we introduce classes of such "blow-up data", i.e., initial functions generating 
finite-time blow-up of solutions. Actually, studying such crucial data will eventually re- 
quire to perform a detailed study of the corresponding elliptic systems with non-Lipschitz 
nonlinearities. 

Parabolic equations. To this end, again beginning with the transformed parabolic equation 
fl2.4p . we consider the separate variable solutions 

(2.29) v{x,t) = {T-t)-^F{x). 

Then F{x) solves the elliptic equation (11.71) in Q, i.e., 

'^F = in 



(2.30) 




on dQ. 



Let F{x) ^ be a solution of the problem (12.301) . which is a key object in the present 
paper. Hence, it follows from (I2.29P that initial data 

(2.31) Vo{x) = CF{x), 

where C 7^ is an arbitrary constant to be scaled out, generate blow-up of the solution 
of ([23D according to (12:291) . 

Hyperbolic equations. Similarly, for the hyperbolic counterpart of (12. 4p . 

(2.32) {ij{v))u = {-ir^'A'^v + v, 
we take initial data in the form 

(2.33) v{x,0) = cF{x) and Vt{x,0) = ciF{x), 

with some constants c and ci such that cci > 0. Then the solution blows up in finite 
time. In particular, choosing 

OO and Ci = ^^S^c^-^ 

with P = — ^^^^^ and B = [^^^j^] " , we have the blow-up solution of (I2.32p in the 
separable form 



v[x 



,t) = {T ~t)^BF{x), where T=(f)''. 



2.3. Blow-up rescaled equation as a gradient system: towards the generic blow- 
up behaviour for parabolic PDEs. Let us briefly discuss another important issue 
associated with the scaling (I2.29p . Consider a general solution v{x,t) of the IBVP for 
( 12. 4p . which blows up first time at t = T. Introducing the rescaled variables 



(2.34) v{x,t) = {T -t)- — w{x,T), T = - ln(T - t) +00 as t 
one can see that w{x, r) then solves the following rescaled equation: 

(2.35) (^(w^))r = (-l)™+^A™u; + w-\ \w\-^w, 
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where on the right-hand side we observe the same operator with a non-Lipschitz nonhn- 
earity as in fll.7l) or fl2.30p . By an analogous argument, fl2.35p is a gradient system and 
admits a Lyapunov functions that is strictly monotone on non-equilibrium orbits: 

(2.36) £{-U\DM' + U^'-^I\M^^) = T^^I kr^k.P > 0. 

Therefore, the corresponding to (12.221) conclusion holds, i.e., all bounded orbits can ap- 
proach stationary solutions only: 

(2.37) uj{wo) C S = {F e H^Q^'in) : (-l)™+iA'"F + F - ^ jFp^F = O}. 

Moreover, since under natural smoothness parabolic properties, uj{wo) is connected and 
invariant [21], the omega-limit set reduces to a single equilibrium provided that S is 
disjoint, i.e., consists of isolated points. Here, the structure of the stationary rescaled set 
S becomes key for understanding blow-up behaviour of general solutions of the higher- 
order parabolic flow (11. ip . 

Thus, the above analysis shows again that the "stationary" elliptic problems (11.70 and 
(I2.30p are crucial for revealing various local and global evolution properties of all three 
classes of PDEs involved. We begin this study with an application of classic variational 
techniques. 

3. Problem "Existence" : VARIATIONAL APPROACH AND COUNTABLE FAMILIES OF 
SOLUTIONS BY LUSTERNIK-SCHNIREL'MAN CATEGORY AND FIBERING THEORY 

3.1. Variational setting and compactly supported solutions. Thus, we need to 
study, in a general multi-dimensional geometry, existence and multiplicity of compactly 
supported solutions of the elliptic problem in (II. 7p . 

Since all the operators in (11.70 are potential, the problem admits the variational setting 
in L^, so the solutions can be obtained as critical points of the following functional: 

(3.1) i?(F) = -i/|D™Fp + i/F2-i/|F|/^, with /?=^G (1,2), 



where, as above, = for even m and D"^ = VA^^ for odd m. In general, we have 
to look for critical points in VF^(]R^)nL^(R^)nL^(M^). Bearing in mind compactly sup- 
ported solutions, we choose a sufficiently large radius i? > of the ball Br and consider 
the variational problem for (13. ip in W^qIBji), where we assume Dirichlet boundary con- 
ditions on Sr = BBr. Then both spaces LF'{Br) and W^^^Br) are compactly embedded 
into W^q{Br) in the subcritical Sobolev range 

(3.2) l<P<Ps = ^ {fi<Ps)- 
In general, we have to use the following preliminary result: 

Proposition 3.1. Let F be a continuous weak solution of the equation ?P such that 

(3.3) F{y) ^0 as \y\ oo. 
Then F is compactly supported in . 
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Notice that continuity of F is guaranteed for < 2m directly by Sobolev embedding 
H"^(M.^) C C(]R^), and, in the whole range (13.21) . by local elliptic regularity theory; see 
necessary embeddings of functional spaces in Maz'ja [271 Ch. 1]. 

Proof. We consider the corresponding parabolic equation with the same elliptic operator, 

— ^ 

(3.4) wt = (-l)'"+^A™u) + w-\w\ "+iu) in x M+, 
with initial data F{y). Setting w = e^w yields the equation 

wt = {-l)'^+^A'^w -e~^^\w\P-^w, where p=^G(0,l), 

where the operator is monotone in L^(M^). Therefore, the Cauchy problem (CP) with 
initial data F has a unique weak solution, [26, Ch. 2]. Thus, (13.41) has the unique so- 
lution w{y,t) = F{y), which then must be compactly supported for arbitrarily small 
t > 0. Indeed, such nonstationary instant compactification phenomena for quasilinear 
absorption-diffusion equations with singular absorption —\u\p~^^u, with p < I, have been 
known since the 1970s and are also called the instantaneous shrinking of the support 
of solutions. These phenomena have been proved for quasilinear higher-order parabolic 
equations with non-Lipschitz absorption terms; see ^9]. □ 

Thus, to revealing compactly supported patterns F{y), we can pose the problem in 
bounded balls that are sufficiently large. Indeed, one can see from (13.11) that, in small 
domains, nontrivial solutions are impossible. 

3.2. L— S theory and direct application of fibering method. The functional (13.11) 
is and is uniformly differentiable and weakly continuous, so we can apply classic 
Lusternik-Schnirel'man (L-S) theory of calculus of variations [221 § 57] in the form of 
the fibering method [SDl EI] , as a convenient generalization of previous versions (TJ of 
variational approaches. 

Namely, following L-S theory and the fibering approach [31], the number of critical 
points of the functional (13. ip depends on the category (or genus) of the functional subset 
on which fibering is taking place. Critical points of E{F) are obtained by spherical fibering 

(3.5) F = r{v)v (r>0), 

where r{v) is a scalar functional, and v belongs to a subset in W^q{Bii) given as follows: 

(3.6) no = {veWl,iBn): Hoiv) ^ - J \D-^v\' + J = l} . 
The new functional 

(3.7) H{r,v) = I r^ - ^rl" J Ivl" 
has the absolute minimum point, where 

(3.8) H'^ = r - rl^-^ J \v\^ = =^ roiv) = {J \vf)^ . 
We then obtain the following functional: 

(3.9) H{v) = H{r,{v),v) = -^rg(^) = (/ kl'^)^- 
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The critical points of the functional fl3.9p on the set (I3.6p coincide with those for 



(3.10) H{v) = Jlvlf^, 

so we arrive at even, non-negative, convex, and uniformly differentiable functional, to 
which L-S theory applies, [251 § 57]; see also [H p. 353]. Following [3T], searching for 
critical points of H on the set Hq one needs to estimate the category p of the set TCq. 
The details on this notation and basic results can be found in Berger [3l p. 378]. Notice 
that the Morse index q of the quadratic form Q in Theorem 6.7.9 therein is precisely the 
dimension of the space where the corresponding form is negatively definite. This includes 
all the multiplicities of eigenfunctions involved in the corresponding subspace. See also 
genus and cogenus definitions and applications to variational problems in [1] and [2] . 
It follows that, by this variational construction, F is an eigenfunction satisfying 

{-l)"^+^A'^F + F - fi\F\ "+iF = 0, 

where yU > is Lagrange's multiplier. Then scaling F i— > yields the original 

equation in (11. 7p . 

For further discussion of geometric shapes of patterns, it is convenient to recall that 
utilizing Berger's version [3l p. 368] of this minimax analysis of L-S category theory [251 
p. 387], the critical values {ck} and the corresponding critical points {vk} are given by 

(3.11) Ck = mi:F(zMk sup^g^ H{v), 

where C Hq are closed sets, and J^k denotes the set of all subsets of the form BS'^^^ C 
TCo, where S"^"^ is a suitable sufficiently smooth (/c — l)-dimensional manifold (say, sphere) 
in Tio and B is an odd continuous map. Then each member of A^^ is of genus at least k 
(available in Tio)- It is also important to remind that the definition of genus [25, p. 385] 
assumes that p{J-') = 1, if no component of JF U JF*, where J-'* = {v : v* = —v E JF}, 
is the reflection of JF relative to 0, contains a pair of antipodal points v and v* = —v. 
Furthermore, /5(JF) = n if each compact subset of JF can be covered by, minimum, n sets 
of genus one. According to (13. lip . 

Ci < C2 < ... < Qo, 

where /q = ^o(-R) is the category of Hq (see an estimate below) such that 

(3.12) /o(^)-^+oo as R^oo. 

Roughly speaking, since the dimension of the sets involved in the construction of A^^ 
increases with k, this guarantees that the critical points delivering critical values (13.110 
are all different. It follows from (13.60 that the category Iq = p(TCo) of the set TCo is equal 
to the number (with multiplicities) of the eigenvalues < 1, 

(3.13) lo = piHo) = ^{Xk < 1} 
of the linear poly-harmonic operator (—I)*" A™ > 0, 

(3.14) (-l)™A™V'fc = hi^k, ipk e Wl.iBR); 
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see [21 p. 368]. Since the dependence of the spectrum on R is, obviously, 
(3.15) Afc(i?) = i?-2™Afc(l), fc = 0,l,2,..., 

we have that the category piTio) can be arbitrarily large for i? 1, and (13.121) holds. We 
fix this in the following: 



Proposition 3.2. The elliptic problem (j i. 7p has at least a countable set of different 



solutions denoted by {F;, I > 0}, each one Fi obtained as a critical point of the functional 
{\3.1\) in W^q{Br) with a sufficiently large R = R{1) > 0. 



Indeed, in view of Proposition 13.11 we choose R^ 1 such that supp Fi C Br. 

3.3. On a model supplying with explicit description of the L— S sequence. As 
we will see shortly, detecting the L-S sequence of critical values for the original functional 
(13. ip is a hard problem, where numerical estimates of the functional will be key. 

However, there exist similar models, for which this can be done much easier. We now 
perform a slight modification in (13. ip and consider the functional 

(3.16) Ei(F) = -i/|£)-Fp + i/F2-i(/F2)t (/? = n±| g (1, 2)) . 
This corresponds to the following non-local elliptic problem: 

(3.17) -(-A)"^F + F-F(/F2)2"^ = (in S^, etc.) 

Denoting by {A^} the spectrum in (I3.14p and by {ipk} the corresponding eigenfunction 
set, we can solve the problem (13.171) explicitly: looking for solutions 

(3.18) F = E(fc>i)Cfc^fc =^ c4-Afc + l- (Eo>i)S')^"'] =0, k = l,2,.... 

The algebraic system in (13.180 is easy and yields precisely the number (I3.13P of various 
nontrivial basic solutions Fi having the form 

(3.19) Fi{y) = ciMy), where |c/-2 = -A^ + 1 > 0, / = 1, 2, Zq- 

3.4. Preliminary analysis of geometric shapes of patterns. The forthcoming dis- 
cussions and conclusions should be understood in conjunction with the results obtained in 
Section [5] by numerical and other analytic and formal methods. In particular, we use here 
the concepts of the index and Sturm classification of various basic and other patterns. 

Thus, we now discuss key questions of the spatial structure of patterns constructed by 
the L-S method. Namely, we would like to know how the genus k of subsets involved in 
the minimax procedure (13.111) can be attributed to the "geometry" of the critical point 
Vk{y) obtained in such a manner. In this discussion, we assume to explain how to merge 
the L-S genus variational aspects with the actual practical structure of "essential zeros 
and extrema" of basic patterns {Fi}. Recall that in the second-order case m = 1, = 1 
this is easy: by Sturm's Theorem, the genus /, which can be formally "attributed" to 
the function Fi, is equal to the number of zeros (sign changes) / — 1 or the number / of 
isolated local extrema points. Though, even for m = 1, this is not that univalent: there 
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are other structures that do not obey the Sturmian order (think about the solution via 
gluing {Fq, Fq, Fq} without sign changes); see more comments below. 

For m > 2, this question is more difficult, and seems does not admit a clear rigorous 
treatment. Nevertheless, we will try to clarify some its aspects. 

Given a solution F of (11. 7p (a critical point of (13.11) ). let us calculate the corresponding 
critical value cp of (I3.10p on the set (13.61) . by taking 

v = CF eHo =^ C = , , ^ 

SO that CF = H{v) = (_j|gj;pVjj.2)^/2 [P = ^)- 
This formula is important in what follows. 

Genus one . As usual in many variational elliptic problems, the ffist pattern Fq (typically, 
called a ground state) is always of the simplest geometric shape, is radially symmetric, 
and is a localized profile such as those in Figure [9l Indeed, this simple shape with a single 
dominant maximum is associated with the variational formulation for Fq: 

(3.21) Fq = r{vQ)vQ, with vq : M H{v) = inf / v G Hq. 

This is (13. lip with the simplest choice of closed sets as points, JF = {v}. 

Let us illustrate why a localized pattern like Fq delivers the minimum to H in (I3.2ip . 
Take e.g. a two-hump structure, 

v{y) = C [vQ{y) +VQ{y + a)\, Ce R, 

with sufficiently large \a\ > diamsuppFo, so that supports of these two functions do not 
overlap. Then, evidently, v G Hq implies that C = and, since (3 G (1,2), 

H{v)=2"-^H{vq)>H{vq). 

By a similar reason, -Fo(y) and VQ{y) cannot have "strong nonlinear oscillations" (see next 
sections for related concepts developed in this direction), i.e., the positive part (Fo)+ must 
be dominant, so that the negative part (-Fq)- cannot be considered as a separate dominant 
1-hump structure. Otherwise, deleting it will diminish H{v) as above. In other words, 
essentially non-monotone patterns such as in Figure [11] or [12] cannot correspond to the 
variational problem (13.211) . i.e., the genus of the functional sets involved is p = 1. 

Radial symmetry of vq is often standard in elliptic theory, though is not straightforward 
at all in view of the lack of the Maximum Principle and moving plane/spheres tools based 
on Aleksandrov's Reflection Principle. We just note that small non-radial deformations 
of this structure, vq i— > vq will more essentially affect (increase) the ffist differential term 
J \D"^vq\'^ rather than the second one in the formula for C in (I3.20p . Therefore, a standard 
scaling to keep this function in Hq would mean taking Cvq with a constant C > 1. Hence, 

HiCvQ) = C^H{vq) ^ C^H{vq) > H{vq), 



so infinitesimal non-radial perturbations do not provide us with critical points of (13.2 ip . 
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For = 1, this shows that Ci cannot be attained at another "positively dominant" 
pattern F+4, with a shape shown in Figure [TST a). See Table 1 below, where for n = 1, 

cf+^ = 1.9488... > C2 = Cf, = 1.8855... > ci = c^o = 1.6203... . 

Genus two . Let now again for simplicity = 1, and let Fq obtained above for the genus 
p = 1 be a simple compactly supported pattern as in Figure [9l By vo{y) we denote the 
corresponding critical point given by f l3.2ip . We now take the function corresponding to 
the difference (15.71) . 

(3.22) V2iy) = vo{y - yo) + v{y + yo)] G Ho (supp vq = [-yo, yo]) , 
which approximates the basic profile Fi given in Figure [TTJ One can see that 

(3.23) H{v2) = 2^H{vo) = 2^ci, 
so that, by fl3llD with k = 2, 

(3.24) ci < C2 < 2^ci. 

On the other hand, the sum as in (15. 6p (cf. Figure [T2l). 

(3.25) V2{y) = ^ [vo{y - yo) + v{y + y^)] e Ho 

also delivers the same value (I3.23P to the functional H. 

It is easy to see that these patterns Fi and F+2,2,+2 as well as -F+4 and many others with 
two dominant extrema can be embedded into a ID subset of genus two on Hq. We show 
such a schematic picture if Figure [H Arrows there indicate the directions of deformations 
of patterns on TYq that can lead to any other profile from such a family, 
embedded into a ID subset of genus two on Ho- 
lt seems that, with such a huge, at least, countable variety of similar patterns, we firstly 
distinguish the profile that delivers the critical value C2 given by (13. lip by comparing the 
values (13.201) for each pattern. The results are presented in Table 1 for n = 1, for which 
the critical values (I3.20p are 

(3.26) ^P = H{CF) = j-jjU^0^j^ (/3 = |). 

The corresponding profiles are shown in Figure El Calculations have been performed with 
the enhanced values Tols= e = 10~^. Comparing the critical values in Table 1, we thus 
arrive at the following conclusion based on this analytical-numerical evidence: 

(3.27) for genus k = 2, the S-L critical value C2 = 1.855... is delivered by Fi. 

Notice that the critical values cp for Fi and -F+2,2,+2 are close by just two percent. 

Thus, by Table 1 the second critical value C2 is achieved at the 1-dipole solution Fi{y) 
having the transversal zero at y = 0, i.e., without any part of the oscillatory tail for y 0. 
Therefore, the neighbouring profile -F-2,3,+2 (see the dotted line in Figure [2]) which has 
a small remnant of the oscillatory tail (see details in Section H]) with just 3 extra zeros, 
delivers another, worse value 

Cf = 1.9268... for F = F_2,3,+2. 
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Figure 1. Profiles Fi = -F-2,1,+2; -^+2,2,+27 etc., can be embedded into a 
ID subset in Hq of genus two. 

Table 1. Critical values of H{v); genus two 



F 


Cf 




Fo 


1.6203.. 


• = Ci 


Fi 


1.8855.. 


• = C2 


F+2,2,+2 


1.9255.. 




F-2,3,+2 


1.9268.. 




F+2,4,+2 


1.9269.. 




F+2,oo,+2 


1.9269.. 




F+A 


1.9488.. 





In addition, the lines from second to fifth in Table 1 clearly show how H increases with 
the number of zeros in between the ±Fo-structures involved. 

Remark: even for m = 1 profiles are not variationally recognizable. Recall that 
for m = 1, i.e., for the ODE (15.51) . the -Fo(?/) profile is not oscillatory at the interface, so 
that the future rule fl3.3ip fails. This does not explain the difference between Fi{y) and, 
say, F+2,0,+2, which, obviously deliver the same critical S-L values by (13. lip . This is the 
case where we should conventionally attribute the S-L critical point to Fi. Of course, for 
m = 1, existence of profiles Fi{y) with precisely / zeros (sign changes) and / + 1 extrema 
follows from Sturm's Theorem. 
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m=2, n=1 : profiles associated with genus p=2 




y 

Figure 2. Seven patterns F{y) indicated in Table 1; m = 2 and n = 1. 

Checking the accuracy of numerics and using fl3.23p . we take the critical values in 
the first and the fifth lines in Table 1 to get for the profile -F+2,00,+25 consisting of two 
independent -Fq's, to within 10~^, 

cf = 2^H{vo) = 2-ici = 1.1892... x 1.6203... = 1.9269... . 

Genus three . Similarly, for A; = 3 (genus p = 3), there are also several patterns that can 
pretend to deliver the L-S critical value C3. These are shown in Figure [31 

The corresponding critical values (13.261) for ?7, = 1 are shown in Table 2, which allows 
to conclude as follows: 

(3.28) k = 3: the S-L critical value C3 = 2.0710... is again delivered by the basic F2. 

All critical values in Table 2 are very close to each other. Again checking the accuracy of 
numerics and taking the critical values Ci in Table 1 and cp for -F+2,00,+2,00,+2 in Table 2, 
consisting of three independent Fq's, yields, to within 10~^, 

Cp = 3^H{vo) = 33ci = 1.31607... x 1.6203... = 2.1324... . 

Note that the S-L category-genus construction (13. lip itself guarantees that all solutions 
{vk} as critical points will be (geometrically) distinct; see [23, P- 381]. Here we stress upon 
two important conclusions: 

(I) First, what is key for us, is that closed subsets in Hq of functions of the sum type in 
(I3.25P do not deliver S-L critical values in (13.110 : 

(II) On the other hand, patterns of the {Fq, Fo}-interaction, i.e., those associated with 
the sum structure as in (I3.25P do exist; see Figures [12] for m = 2; and 
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Table 2. Critical values of H{v); genus three 

F Cf 

F2 2.0710... = C3 

-f+2,2,+2,2,+2 2.1305... 
-^+2,3,-2,3,4-2 2.1322... 
-^+2,00, +2,00, -f2 2.1324... 
F+6 2.1647... 



(Ill) Hence, these patterns (different from the basic ones {Fi}) as well as many others 
are not obtainable by a direct S-L approach. Therefore, we will need another version of 
S-L and fibering theory, with different type of decomposition of functional spaces to be 
introduced below. 

Genus k . Similarly taking a proper sum of shifted and reflected functions ±Vo{y ± lyo), 
we obtain from (13. lip that 

(3.29) Ck-i <Ck< k^ci. 

Conclusions: conjecture and an open problem. As a conclusion, we mention that, regard- 
less such close critical values cp, the above numerics confirm that there is a geometric- 
algebraic way to distinguish the S-L patterns delivering (13.111) . It can be seen from (13.261) 
that, destroying internal oscillatory "tail" , or even any two-three zeros between two Fo-like 
patterns in the complicated pattern F{y), 

(3.30) decreases two main terms — J and J |F|t in cf in (13.201) . 
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Recall that precisely these terms in the ODE 

F^"^) = -\F\-^F + ... (see dED forn = 1) 

are responsible for formation of the tail as shown in Section HI while the F-term, giving 
J F"^, is negligible in the tail. Decreasing both terms, i.e., eliminating the tail in between 
-Fq's, will decrease the value cp, since in (13.201) the numerator gets less and the denominator 
larger. Therefore, composing a complicated pattern Fi{y) from several elementary profiles 
looking like Fo{y), by using {k — 1) -dimensional manifolds of genus k, we follow 

Formal Rule of Composition (FRC) of patterns: performing maximization of H{v) 
of any [k — 1)- dimensional manifold T G M.^, 

(3.31) the S-L point Fk~i{y) is obtained by minimizing all internal tails and zeros, 

i.e., making the minimal number of internal transversal zeros between single structures. 

Regardless such a simple variational-oscillatory meaning (13.301) of this FRC, we do not 
know how to give to such a rule a rigorous sounding. 

Concerning the actual critical S-L points, we end up with the following conjecture, 
which well- corresponds to the FRC (I3.3ip : 

Conjecture [31.4. For N = 1 and any m > 2, the critical S-L value k > 1, 

is delivered by the basic pattern Fk^i, obtained by minimization on the corresponding 
[k — 1)- dimensional manifold G M.k, which is the interaction 

(3.32) Fk-i = (-l)'^"M+Fo, -Fo, +Fo, {-if-^Fo}, 

where each neighbouring pair {Fq, — Fq} or {— Fo,Fo} has a single transversal zero in 
between the structures. 

We also would like to formulate the following conclusion that is again associated with 
the specific structure of the L-S construction (13.111) over suitable subsets as smooth 
{k — l)-dimensional manifolds of genus k: 

Open problem [31.5. For N = 1 and m > 2, there are no any purely geometric-topology 
arguments establishing that the conclusion of Conjecture [^4 holds. Naturally, the same 
remains true in M^. 

In other words, we claim that the metric "tail" -analysis of the functionals involved in 
the FRC (I3.3ip cannot be dispensed with by any geometric-type arguments. Actually, 
the geometric analysis is nicely applied for m = 1 and this is perfectly covered by Sturm's 
Theorem on zeros for second-order ODEs. If such a Sturm Theorem is nonexistent, this 
emphasizes the end of geometry-topology (or purely homotopy if tails are oscillatory) 
nature of the variational problem under consideration. 

On patterns in R^ . In the elliptic setting in M^, such a clear picture of basic patterns 
Fi{y) obtained via (13.111) with A; = / -|- 1 is not available. As usual in elliptic theory, nodal 
structure of solutions in is very difficult to reveal. Nevertheless, we strongly believe 
that the L-S minimax approach (13. lip also can be used to detect the geometric shape of 
patterns in the A^-dimensional geometry. 
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For instance, it is most plausible that the 1-dipole profile -Fi(y), which is not radially 
symmetric is essentially composed from two radial FQ{y)-tjpe profiles via an {—Fq,Fq}- 
interaction {q.v. Figure [15] in ID). Therefore, Fi{y) has two dominant extrema, in a 
natural way, similar to the second eigenfunction 

(3.33) d^e-l^'I'/^-^/ie-IJ'l'/^ in 
of the self-adjoint second-order Hermite operator 

H2 = A + lyV + fL 

Such a comparison assumes a bifurcation phenomenon at n = from eigenfunctions of the 
linearized operator; see applications in [T51 Appendix] to second-order quasilinear porous 
medium operators. Compactification of the pattern fl3.33p and making it oscillatory at 
the interface surface would lead to correct understanding what the 1-dipole profile looks 
like, at least, for small n > 0. Similar analogy is developed for all odd patterns F2k+i- 
For instance, F^ly) has a dominant "topology" similar to the fourth eigenfunction of H2, 

which has precisely four extrema on the j/i-axis, etc. 

Concerning even patterns F2k, we believe that the above L-S algorithm leads to simple 
radially-symmetric solutions of (11.71) . i.e., solutions of ODEs; see further comments below. 

Remark: on radial and ID geometry. Of course, the elliptic equation (11.71) also 
admits a countable family of radially symmetric solutions 

{FrWvD, I = 0,2,4,.:} 

satisfying the corresponding ODE. These are constructed in a similar manner by L-S and 
fibering theory. In ID, this gives the basic set / = 0, 1, 2, ...} that was described in 

the previous sections. We expect that the first members of all three families coincide, 

(3.34) Fo = F^"^ = F^^. 

Further correspondence of the L-S spectrum of patterns in Proposition 13.21 and the ID 
one / > 0} is discussed in [181 § 4]. 

4. Problem "Oscillations" : local oscillatory structure of solutions 

CLOSE TO INTERFACES AND PERIODIC CONNECTIONS WITH SINGULARITIES 

As we have seen, the first principal feature of the ODEs (II. 8p (and the elliptic coun- 
terparts) is that these admit compactly supported solutions. Indeed, all interesting for 
us patterns have finite interfaces. This has been proved in Proposition 13.11 in a general 
elliptic setting. 

Therefore, we are going to study typical local behaviour of the solutions of (II. Sp close to 
the singular points, i.e., to finite interfaces. We will reveal extremely oscillatory structure 
of such behaviour to be compared with global oscillatory behaviour obtained above by 
variational techniques. 
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The phenomenon of oscillatory changing sign behaviour of solutions of the Cauchy 
problem has been detected for various classes of evolution PDEs; see a general view in 
[20| Ch. 3-5] and various results for different PDEs in [HI |12l |13] . For the present 2mth- 
order equations, the oscillatory behaviour exhibits special features to be revealed. We 
expect that the presented oscillation analysis makes sense for more general solutions of 
the parabolic equation (11. ip and explains their generic behaviour close to the moving 
interfaces. 

4.1. Autonomous ODEs for oscillatory components. Assume that the finite inter- 
face of F{y) is situated at the origin y = 0, so that we can use the trivial extension 
F{y) = for y < 0. We then are interested in describing the behaviour of solutions as 
y — i> O"*", so we consider the ODE fll.Sp written in the form 

(4.1) = (_l)m^ ^ > jT^^O) = (a=^G(0,l)). 

In view of the scaling structure of the first two terms, for convenience, we perform extra 
rescaling and introduce the oscillatory component 99(5) of F by 

(4.2) F{y) = y^^{s), where s = \ny and 7 = ^ = ^^^i^^. 

Therefore, since s —00 (the new interface position) as ?/ — > 0^, the monotone function 
y'^ in (14. 2 p plays the role of the envelope to the oscillatory function F{y). Substituting 
(14. 2 p into (14.11) yields the following equation for ip: 

(4.3) P2M = (-1)"^+Vr> + (-l)'"e2™V- 

Here {Pk, k > 0} are linear differential operators defined by the recursion 

(4.4) Fk+ii^p) = {Pkiv))' + (7 - k)Pk{^) for = 0, 1, ... , Po(v^) = ^- 
Let us present the first five operators, which are sufficient for further use: 

Fi((/j) = </?' + 7¥?; 

P2(¥^) = <^" + (27-l)v^' + 7(7-l)<^; 

P^{ip) = + 3(7 - \y + (37' - 67 + iy + 7(7 - 1)(7 - 2)</^; 
P^{ip) = <^(^) + 2(27 - 3)ip"' + (67' - I87 + 11)^" 

+ 2(27^ - 97" + II7 - 3)ip' + 7(7 - 1)(7 - 2)(7 - 3)ip; 
P^iip) = (^(^) + 5(7 - 2)(^(^) + 5(27^ - 87 + 7)ip"' + 5(7 - 2) (272 - 87 + 5)ip" 

+ (57^ - 407^ + 1057^ - IOO7 + 24)(^' + 7(7 - 1)(7 - 2) (7 - 3) (7 - 4)(^; etc. 

According to (14.21) . the interface at y = now corresponds to s = —00, so that (14. 3 p is 
an exponentially (as s — *• —00) perturbed autonomous ODE 

(4.5) P2„^(¥^) = (-l)^'"+%|-> in M (a = ^), 

which we will concentrate upon. By classic ODE theory [8], one can expect that for 
s <^ — 1 typical (generic) solutions of (14. 3 p and (14. 5 p asymptotically differ by exponentially 
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small factors. Of course, we must admit that (14. 5 p is a singular ODE with a non-Lipschitz 
term, so the results on continuous dependence need extra justification in general. 
Thus, in two principal cases, the ODEs for the oscillatory component ip{s) are 

(4.6) m = 2: P,{if) = ~\if\-"ip, 

(4.7) m = 3: P^^icp) = +\ip\-''ip, 

that exhibit rather different properties because comprise even and odd m's. For instance, 
(14. 5 p for any odd m > 1 (including (14.71) ) has two constant equilibria, since 

7(7- 1)... (7- (2m- !))(/?= Iv^r^'V' =^ 
^'^ y?(s) = ±(/?o = ±[7(7- l)---(7- (2m- 1))]-^ for all n > 0. 

For even m including (14.60 . such equilibria for (14. 5 P do not exist at least for n G (0, 1]. 
We will show how this affects the oscillatory properties of solutions for odd and even m's. 

4.2. Periodic oscillatory components. We now look for periodic solutions of (14. Sp . 
which are the simplest nontrivial bounded solutions that can be continued up to the 
interface at s = —00. Periodic solutions together with their stable manifolds are simple 
connections with the interface, as a singular point of the ODE (II. 8p . 

Note that (14. 5 p does not admit variational setting, so we cannot apply well developed 
potential [28l Ch. 8] (see a large amount of related existence-nonexistence results and 
further references therein), or degree [211 [25] theory. For m = 2, the proof of existence 
of can be done by shooting; see [HI § 7.1] that can be extended to m = 3 as well. 
Nevertheless, uniqueness of a periodic orbit is still open, so we conjecture the following 
result supported by various numerical and analytical evidence (cf. f2U\ § 3.7]): 

Conjecture [Hi. For any m > 2 and a G (0, 1], the ODE ( |^.5p admits a unique nontrivial 
periodic solution (p^:{s) of changing sign. 

4.3. Numerical construction of periodic orbits for m = 2. Numerical results clearly 
suggest that (14. 6 p possesses a unique periodic solution f*{s), which is stable in the direc- 
tion opposite to the interface, i.e., as s ^ +00; see Figure [H The proof of exponential 
stability and hyperbolicity of is straightforward by estimating the eigenvalues of the 
linearized operator. This agrees with the obviously correct similar result for n = 0, 
namely, for the linear equation (14. ip for a = 

(4.9) F(^) = -F as 2/ -> -00. 

Here the interface is infinite, so its position corresponds to y = —00. Indeed, setting 
F{y) = gives the characteristic equation and a unique exponentially decaying pattern 

(4.10) /i^ = -l =^ F(?/) ~e^ [Acos(^) +Ssin(^)] as 2/-^ -00. 

Continuous dependence on r?, > of typical solutions of (14.50 with "transversal" zeros only 
will continue to be key in our analysis, that actually means existence of a "homotopic" 
connection between the nonlinear and the linear (n = 0) equations. The passage to 
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(a) n = 2 (b) n = 4 

Figure 4. Convergence to the stable periodic solution of ()4.6p for n = 2 (a) 
and n = 4 (b). 

the limit n — >^ in similar degenerate ODEs from thin film equations (TFEs) theory is 
discussed in [HI § 7.6]. 

The oscillation amplitude becomes very small for n ~ 0, so we perform extra scaling. 

Limit n — » . This scaling is 

(4.11) ^{s)={^)-^v), where = f , 
where $ solves a simpler limit binomial ODE, 

TL 

(4.12) e-''(e''$)(^) = + 4$'" + 6$" + 4$' + $ = -|$|"^$. 

The stable oscillatory patterns of fl4.12p are shown in Figure O For such small n in Figure 
O^a) and (b), by scaling fl4.1ip . the periodic components (p^ get really small, e.g., 

max|(^,(s)| ~ 3 ■ lO-^(f)" ~ 3 ■ 10-30 for n = 0.2 in (a), 
and max |y9=K(s)| ~ 10^^^ for n = 0.08 in (b). 

Limit n — » Then a — >■ 1, so the original ODE ( 14. GP approaches the following equation 
with a discontinuous sign-nonlinearity: 

(4.13) + lOipZ + 35<^'4 + 50^'^ + 24(^00 = -sign 
This also admits a stable periodic solution, as shown in Figure [61 

4.4. Numerical construction of periodic orbits for m = 3. Consider now equation 
(14. 7p that admits constant equilibria (14.80 existing for all n > 0. It is easy to check that 
the equilibria ±ipo are asymptotically stable as s — +oo. Then the necessary periodic 
orbit is situated in between of these stable equihbria, so it is unstable as s ^ +oo. 

Such an unstable periodic solution of (14.70 is shows in Figure [7] for n = 15, which is 
obtained by shooting from s = with prescribed Cauchy data. 
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(a) n = 0.2 (b) n = 0.08 

Figure 5. Stable periodic oscillations in the ODE (I4.12P for n = 0.2 (a) and 
n = 0.08 (b). 
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Figure 6. Convergence to the stable periodic solution of (14.130 (n = +00). 

As for m = 2, in order to reveal periodic oscillations for smaller n (actually, there is a 
numerical difficulty already for n < 4), we apply the scaling 

6. 

(4.14) v9(s) = (1) "$(77), where V = ^- 

This gives in the limit a simplified ODE with the binomial linear operator, 

(4.15) e-''(e''<l>)(^) = <^^^^ + 6<l>^^^ + 15<l>(^) + 20$'" + 15$" + 6$' + $ = 

Figure [8] shows the trace of the periodic behaviour for equation (14. 151) with n = i. Ac- 
cording to scaling (14.141) . the periodic oscillatory component v?*(s) gets very small, 

max \ip^\ ~ 1.1 X 10"^*^ for n = 0.5. 
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Figure 7. Unstable periodic behavior of the ODE (fiTTl) for n = 15. Cauchy 
data are given by (^(0) = 10"^, ip'iO) = ip"'{0) = ... = ip'^^\0) = 0, and ip"iO) = 
-5.0680839826093907... x lO""^ . 
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Figure 8. Unstable periodic behavior of the ODE (|4.15[) for n = ^. Cauchy 
data are given by (^(0) = lO""^, ip'{0) = ip"'{0) = ... = ip^^\0) = 0, and ip"{0) = 
-9.456770333415... x lO""^ . 

A more detailed study of the behaviour of the oscillatory component as n — > is available 
in da § 12]. 

The passage to the limit n +oo leads to the equation with discontinuous nonlinearity 
that is easily obtained from (14. 7p . This admits a periodic solution, which is rather close 
to the periodic orbit in Figure [7| obtained for n = 15. 

We claim that the above two cases m = 2 (even) and m = 3 (odd) exhaust all key types 
of periodic behaviours in ODEs like (11.81) . Namely, periodic orbits are stable for even m 
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and are unstable for odd, with typical stable and unstable manifolds as s ^ ±00. So we 
observe a dichotomy relative to all orders 2m of the ODEs under consideration. 



5. Problem "Numerics" : NUMERICAL CONSTRUCTION AND FIRST CLASSIFICATION 
OF BASIC TYPES OF LOCALIZED BLOW-UP OR COMPACTON PATTERNS FOR m = 2 

We need a careful numerical description of various families of solutions of the ODEs 
(11. 8p . In practical computations, we have to use the regularized version of the equations. 



m 



n 

(5.1) (^-i)^F^^"'^ = F - [e"^ + F'^y^^F 

which, for £ > 0, contains smooth analytic nonlinearities. In numerical analysis, we 
typically take e = 10""^ or, at least, 10~^ which is sufficient to revealing global structures. 

It is worth mentioning that detecting in Section H] a highly oscillatory structure of 
solutions close to interfaces makes it impossible to use well-developed homotopy theory 
[22| Hi] that was successfully applied to another class of fourth-order ODEs with coercive 
operators; see also Peletier-Troy [29]. Roughly speaking, our non-smooth problem cannot 
be used in a homotopy classification, since the oscillatory behaviour close to interfaces de- 
stroys such a standard homotopy parameter as the number of rotations on the hodograph 
plane {F, F'}. Indeed, for any solution of (II. 8p . the rotation number about the origin is 
always infinite. Then as F — 0, i.e., as y ^ ±00, the linearized equation is (14. 9p . which 
admits the oscillatory behaviour (14.101) . 

5.1. Fourth-order equation: m = 2. We will describe main families of solutions. 
First basic pattern and structure of zeros. For m = 2, (11.81) reads 

(5.2) F(^) = F -\f\ in M. 

We are looking for compactly supported patterns F (see Proposition 13. ip satisfying 

meassupp F > 2/2=,,, where i?* > ^ is the first 

(5.3) ^ 

positive root of the equation tanh R = — tan R. 

In Figure [9], we show the first basic pattern called the Fo^y) for various n G [j^, 100]. 
Concerning the last profile n = 100, note that (15. 2p admits a natural passage to the limit 
n —>■ +00 that gives the ODE with a discontinuous nonlinearity, 

. X m (F-l for F>0, 

A unique oscillatory solution of (15.40 can be treated by an algebraic approach; cf. [HI 
§ 7.4]. For n = 1000 and n = +00, the profiles are close to that for n = 100 in Figure O 
The profiles in Figure M are constructed by MATLAB with extra accuracy, where e in 
(15. ip and both tolerances in the bvp4c solver have been enhanced and took the values 
e = 10^^ and Tols = 10^^. This allows us also to check the refined local structure of 
multiple zeros at the interfaces. Figure [10] corresponding to r?, = 1 explains how the zero 
structure repeats itself from one zero to another in the usual linear scale. 
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m=2: basic pattern F^{y) for various n 




Figure 9. The first (stable) solution Fo{y) of ()5.2p for various n. 



.3 m=2, n=1: basic pattern F (y), oscillations enlarged 




m=2, n=1: basic pattern F (y), oscillations enlarged 




(a) scale lO"^ (b) gcale lO"'^ 

Figure 10. Enlarged zero structure of the profile -Fo(y) for n = 1 in the linear scale. 



Basic countable family: approximate Sturm 's property. In Figure [TH we show the basic 
family denoted by {F/, I = 0, 1,2, ...} of solutions of (15. 2p for n = 1. Each profile Fi{y) 
has precisely / + 1 "dominant" extrema and / "transversal" zeros; see further discussion 
below and [161 § 4] for other details. It is important that all the internal zeros of Fi{y) 
are clearly transversal (obviously, excluding the oscillatory end points of the support). In 
other words, each profile Fi is approximately obtained by a simple "interaction" (gluing 
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m=2, n=1 : basic family {F|} 

1.5 



F(y)i 



0.5 


-0.5 



-1 
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-20 -15 -10 -5 5 10 15 20 

y 

Figure 11. The first five patterns of the basic family {FJ of the ODE (jO]) 
for n = 1. 

together) of / + 1 copies of the first pattern ±Fo taking with necessary signs; see further 
development below. 

Actually, if we forget for a moment about the complicated oscillatory structure of solu- 
tions near interfaces, where an infinite number of extrema and zeros occur, the dominant 
geometry of profiles in Figure [11] looks like it approximately obeys Sturm's classic zero 
set property, which is true rigorously for m = 1 only, i.e., for the second-order ODE 

n 

(5.5) F" = -F+\F\ "+^F in M. 

For (15. 5p . the basic family {Fi} is constructed by direct gluing together the explicit pat- 
terns (11.51) . i.e., ±Fq. Therefore, each Fi consists of precisely /-|-1 patterns (II. 5p (with signs 
±Fo), so that Sturm's property is clearly true. In Section [31 we presented some analytic 
evidence showing that precisely this basic family {F/} is obtained by direct application of 
L-S category theory. 

5.2. Countable family of {Fq, Fo}-interactions. We now show that the actual non- 
linear interaction of the two first patterns +FQ{y) leads to a new family of profiles. 

In Figure [121 n = 1, we show the first profiles from this family denoted by {F+2,k,+2}, 
where in each function F+2,fc,+2 the multiindex a = {+2, k, +2} means, from left to right, 
+2 intersections with the equilibrium -|-1, then next k intersections with zero, and final 
+2 stands again for 2 intersections with -|-1. Later on, we will use such a multiindex 
notation to classify other patterns obtained. 

In Figure [131 we present the enlarged behaviour of zeros explaining the structure of 
the interior layer of connection of two profiles ~ +Fo{y). In particular, (b) shows that 
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m=2, n=1 : family {F^^ ^ J 

1.5 



F(y) 

0.5 





-15 -10 -5 5 10 15 

y 

Figure 12. First patterns from the family {F+2,k,+2} of the {Fo,Fo}- 
interaction; n = 1. 

there exist two profiles F+2,6,+2, these are given by the dashed line and the previous one, 
both having two zeros on [—1,1]. Therefore, the identification and the classification of 
the profiles just by the successive number of intersections with equilibria and ±1 is not 
always acceptable (in view of a non-homotopic nature of the problem), and some extra 
geometry of curves near intersections should be taken into account. In fact, precisely 
this proves that a standard homotopy classification of patterns is not consistent for such 
non-coercive and oscillatory equations. Anyway, whenever possible without confusion, we 
will continue use such a multiindex classification, though now meaning that in general 
a profile F^. with a given multiindex a may denote actually a class of profiles with the 
given geometric characteristics. Note that the last profile in Figure ITU is indeed F+2,6,+25 
where the last two zeros are seen in the scale ~ 10~^ in Figure dH Observe here a clear 
non-smoothness of two last profiles as a numerical discrete mesh phenomenon, which 
nevertheless does nor spoil at all this differential presentation. 

In view of the oscillatory character of Fo{y) at the interfaces, we expect that the family 
{F+2,fc,+2} is countable, and such functions exist for any even /c = 0, 2, 4, ... . Then k = +00 
corresponds to the non-interacting pair 

(5.6) Fo{y + yo) + Fo{y - yo), where supp Fo{y) = [-yo, Vo]- 

Of course, there exist various triple {Fq, Fq, Fq} and any multiple interactions {Fq, Fq} 
of k single profiles, with different distributions of zeros between any pair of neighbours. 

5.3. Countable family of {— Fq, Foj-interactions. We now describe the interaction 
of —Fo{y) with Fo{y). In Figure [T5l n = 1, we show the first profiles from this family 
denoted by {F_2,a:,+2}, where for the multiindex cr = {— 2, A;, -|-2}, the first number —2 
means 2 intersections with the equilibrium —1, etc. The zero structure close to y = is 
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(a) zeros: scale 10^^ (b) zeros: scale 10"'' 

Figure 13. Enlarged middle zero structure of the profiles Fj^2,k,+2 from Figure [T2l 



fii=2, n=1 : family {F^^^^^}' further enlarged zero structure 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 
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Figure 14. Enlarged middle zero structure of the profiles F+2,6,-1-2 from Figure [T2l 

presented in Figure [IB It follows from (b) that the first two profiles belong to the class 
-F_2,i,2, i-e., both have a single zero for y ^ 0. The last solution shown is F_2,5,-i-2- Again, 
we expect that the family {F_2,a:,+2} is countable, and such functions exist for any odd 
k = 1,3,5, and k = +00 corresponds to the non-interacting pair 

(5.7) ~Fo{y + yo) + Fo{y - yo) (supp Fo{y) = [-yo, yo]). 

There exist families of an arbitrary number of interactions such as {±Fq,±Fq, ...,±Fq} 
consisting of any k >2 members. 

5.4. Periodic solutions in M. Before introducing new types of patterns, we need to de- 
scribe other non-compactly supported solutions in M. As a variational problem, equation 
(15. 2p admits an infinite number of periodic solutions; see e.g. [281 Ch. 8]. In Figure [T7I 
for n = 1, we present a special unstable periodic solution obtained by shooting from the 
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m=2, n=1: family {F } 




Figure 15. First four patterns from the family {F_2^k,+2} of the {— Foj-Foj- 
interactions; n = 1. 



m=2, n=1: family {F ) 



m=2, n=1: family {F ) 



F(y) 




F(y) 



(a) zeros: scale 10 ^ (b) zeros: scale 10 

Figure 16. Enlarged middle zero structure of the profiles F_2^k,+2 from Figure [T5l 



origin with conditions 

F(0) = 1.5, F'(0) = F"'(0) = 0, F"(0) = -0.3787329255... . 
We will show next that precisely the periodic orbit F^:{y) with 
(5.8) F,(0) ^ 1.535... 

plays an important part in the construction of other families of compactly supported 
patterns. Namely, all the variety of solutions of (15. 2p that have oscillations about equilibria 
±1 are close to ±F*(?/) there. 
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Figure 18. Two families of solutions of (j5.2p for n = 1; Fj^2k{y) (a) and 

F+k,l~m,l,+k (b). 



5.5. Family {F+2fc}- Such functions F^2k ^oi k > 1 liave 2k intersection with the single 
equilibrium +1 only and have a clear "almost" periodic structure of oscillations about; 
see Figure [T8](a). The number of intersections denoted by +2k gives an extra Strum index 
to such a pattern. In this notation, F+2 = Fq. 

5.6. More complicated patterns: towards chaotic structures. Using the above 
rather simple families of patterns, we claim that a pattern (possibly, a class of patterns) 
with an arbitrary multiindex of any length 



(5.9) 



a = {±cri, 0-2, ±0-3, (T4, ±ai} 
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m=2, n=1 : pattern F^^ 2+22+6 rr\=2, n=1 : pattern F^^ 




(c) cr = {+2, 2, +4, 2, +2, 1, -4} (d) cr = {+6, 3, -4, 2, -6} 

Figure 19. Various patterns for ([O]) for n = 1. 



can be constructed. Figure [T8l(b) shows several profiles from the family with the index 
a = {+k,l,—m,l,+k}. In Figure [T9| we show further four different patterns, while in 
Figure [20], a single most complicated pattern is presented, for which 

(5.10) a = {-8, 1, +4, 1, -10, 1, +8, 1, 3, -2, 2, -8, 2, 2, -2}. 

All computations are performed for n = 1 as usual. Actually, we claim that the multiindex 
(15. 9p can be arbitrary and takes any finite part of any non-periodic fraction. Actually, 
this means chaotic features of the whole family of solutions {F„}. These chaotic types of 
behaviour are known for other fourth-order ODEs with coercive operators, [291 P- 198]. 
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m=2, n=1: pattern ^ _^^,^,,^,^_,ax,^^^_^.2.-B^^^-2 




6. Problem "Numerics" : PATTERNS IN ID IN HIGHER-ORDER CASES, m > 3 

The main features of the pattern classification by their structure and computed critical 
values for m = 2 in the previous section can be extended to arbitrary m > 3 in the ODEs 
(11. 8p for = 1, so we perform this in less detail. 

In Figure [211 fo^' the purpose of comparison, we show the first basic pattern Fo(y) 
for n = 1 in four main cases m = 1 (the only non-negative profile by the Maximum 
Principle known from the 1970s [38], [371 Ch. 4]), and 2, 3, 4. Next Figure [221 explains 
the oscillatory properties close to the interface. It turns out that, for m = 4, the solutions 
are most oscillatory, so we it is convenient to use this case for illustrations. 

In the log-scale, the zero structure is shown in Figure [23] for m = 2,3, and 4 (n = 1). 
For m = 4 and m = 3, this makes it possible to observe a dozen of oscillations that well 
correspond to the oscillatory component analytic formulae (14. 2 p close to interfaces. For 
the less oscillatory case m = 2, we observe 4 reliable oscillations up to lO"^'', which is our 
best accuracy achieved. 

The basic countable family satisfying approximate Sturm's property has the same topol- 
ogy as for m = 2 in Section [5l and we do not present such numerical illustrations. 

In Figure [23] for m = 3 and n = 1, we show the first profiles from the family {F+2fc}, 
while Figure [23] explains typical structures of F^2,k,+2} for m = 4, n = 1. In Figure [2B] for 
m = 4 and tj, = 1, we show the first profiles from the family {F+2,fe,-2}- 

Finally, in Figure [27] for comparison, we present a complicated pattern for m = 3 and 
4 (the bold line), n = 1, with the index 

(6.1) a = {-8, 3, +4, A;, -10, 1, +8,/, -12}. 

Both numerical experiments were performed starting with the same initial data. As a 
result, we obtain quite similar patterns, with the only difference that, in (16.11) . k = 1, 
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First pattern F =F for nn=1 ,2,3,4 
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Figure 21. The first solution Fo(y) of ([THl), iV = 1, n = 1, for m = 1,2,3,4. 



First pattern F. for m=1 ,2,3,4: oscillations near interface _3 First pattern F. for m=1 ,2,3,4: oscillations enlarged 

" X 10 " 




y y 

(a) scale 10-2 (b) scale IQ-^ 



Figure 22. Enlarged zero structure of the profile Fo{y) for n = 1 from Figure 
[2H the linear scale. 



/ = 3 for m = 3, and for more oscillatory case m = 4, the number of zeros increase, so 
now k = 3 and / = 5. 

Our study of other key aspects of these challenging elliptic problems will be soon con- 
tinued in [18]. 
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m=2,3,4, n=1 : oscillations in log-scale 



log|F(y)| 




10 15 

y 

Figure 23. Behaviour of -Fo(y) for n = 1, for m = 2, 3,4; the log-scale. 

nn=3, n=1: patterns F^^^, k=1 ,2,3,4,5,6,7 




Figure 24. The first seven patterns from the family {F+2fc}; m = 3 and n = 1. 
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m=4, n=1 : seven symmetric patterns F^^ for k=2,4,6 
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